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Abstract 

We extend existing calculations of the electromagnetic dipole operator contribution to 
the total decay rate and the photon energy spectrum of the decay B — > X^^ at O(a^) 
by working out the exact dependence on the charm quark mass. 



1 Introduction 



As a flavor changing neutral current process the inclusive decay B X^'-f is loop-induced 
and therefore highly sensitive to new degrees of freedom beyond the Standard Model. To 
tap the full potential of this decay channel in deriving constraints on the parameter space 
of new physics models both the experiments and the Standard Model calculations should be 
known as accurately as possible. 

On the experimental side, the latest measurements by Belle and BABAR are reported 
in [1,2], and the world average performed by the Heavy Flavor Averaging Group [3] for 
E^, > 1.6 GeV reads 



where the errors are statistical, systematical, due to the extrapolation to the common lower- 
cut in the photon energy, and due to the B Xd'j contamination, respectively. 

In order to compete with the given experimental accuracy the theoretical prediction of the 
B —* branching ratio has to be known at the next-to-next-to-leading order (NNLO) level. 
There have been great efforts of several groups within the last few years to achieve this goal. 
The three- loop dipole operator matching was found in [4], the three-loop mixing of the four- 
quark operators in [5], and the three- loop mixing of the dipole operators was calculated in [6]. 
Furthermore, the four-loop mixing of the four-quark operators into the dipole operators 
was calculated in [7]. The two-loop matrix elements of the electromagnetic dipole operator 
together with the corresponding bremsstrahlung terms can be found in [8-11]. The three-loop 
matrix elements of the four-quark operators were found in [12] within the so-called large- 
(3q approximation, and a calculation that goes beyond this approximation by employing an 
interpolation in the charm quark mass rric was presented in [13]. The combination of all these 
individual contributions culminated in a first estimate of the B Xs'j branching ratio at 
0{al) [14]. For > l.Q GeV it reads 



Here, we should mention that there are several perturbative and non-perturbative effects that 
have not been considered when deriving this estimate. Some of these are already available in 
the literature: the four- loop mixing oi Oi, . . . ,0q into Os [7]; the bremsstrahlung contribu- 
tions of the (O25O2), (02,07) and (O7, 08)-interferences at 0(a^/3o) [15]; those (02,02) and 
(O2, O7) contributions which are due to the renormalization of rric [16] (written as expansions 
in nic/mh); photon energy cut-off related effects [17-21]; and estimates for the 0(«sAqcd/''^6) 
corrections [22]. Other effects are unknown at the moment, like the complete virtual- and 
bremsstrahlung contributions to the (07,03)- and (Og, 08)-interferences at O(tt^) ^; the ex- 
act mc-dependence of various matrix elements beyond the large /3o-approximation, in order 
to improve (or even remove) the uncertainty due to the interpolation in rric [13]; and the emis- 
sion process of photons from four-quark operators at tree-level. The individual contributions 
listed above are all expected to remain within the uncertainty given in (1.2), nevertheless 
they should be taken into account in future updates. 

In the present paper we extend the calculations of the (O7, 07)-interference contribution 
performed in [8-11] to include the charm quark mass at its physical value. Since the results 

-"^In (1.2) the virtual corrections to the (O7, Og) -interference at 0{a'^Po) were taken into account. 



Bt{B X,7) = (3.55 ± 0.24^°;?^ ± 0.03) x 10"^ , 



(1.1) 



Bt{B X,7) = (3.15 ± 0.23) x 10 



,-4 



(1.2) 
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given there were presented for N}j heavy quarks with masses equal to mf, and A'^^ hght quarks 
with masses equal to zero, exact results for the (O7, O7) contribution are only available for 
two extreme cases. Either we can set A'^^ = 2 and A^^^: = 3, that is the charm quark mass 
is equal to m^, or we can set = 1 and Nl = 4, that is the charm quark is considered 
to be massless. The Nh heavy quarks are of course kinematically not allowed to appear in 
the final state, and hence the first possibility corresponds qualitatively to the experimental 
situation since events with charmed hadrons are not included on the experimental side. 
However, since rric ~ fnh/4: in reality, one could choose as well the second possibility, but in 
this case the inclusion of contributions from the cc production is required on the theoretical 
side in order to get rid of infrared divergences. As argued in [10] we expect the true result 
to lie somewhere in between. To check this statement, we calculate the exact charm quark 
mass dependence of the (O7, 07)-interference contribution to the photon energy spectrum 
dT{b XP''''*°°'''7) /dE^ and the total decay width T{b ^ XP'^''*™^'''7), excluding charm quarks 
in the final state. The impact of the exact mc-dependence on the branching ratio will be 
taken into account together with other new contributions in a forthcoming analysis. 

The organization of this paper is as follows. In section 2 we describe briefly the calcula- 
tion of the relevant Feynman diagrams and present our final results for the photon energy 
spectrum and the total decay width. Furthermore, we comment on the numerical importance 
of our result on the branching ratio Br(i? — > Xs'~f) in this section. 



2 Results for the charm quark contribution 

Within the low-energy effective theory the partonic b — Xs'~f decay rate can be written as 

T{b xr'-i)E,>E. = ^>^-^"^[^'^< \v,,v,:\^ ^^'(/^) (^) ^^^(^0, /i) , (2.1) 

where rrtb and mb{fi) denote the pole and the running MS mass of the b quark, respectively, 
Cj^^(yu) the effective Wilson coefficients at the low-energy scale, and Eq the energy cut in the 
photon spectrum. As already anticipated in the introduction, we will focus on the function 
G77{Eq, fi) corresponding to the self-interference of the electromagnetic dipole operator 

^7 = TTT^ (sa^'PRb) F,, . (2.2) 

lovr^ 

More precisely, we extend the calculations performed in [8-11] to include the effects of a 
massive charm quark. To this end we calculate the cuts of the b quark self-energies displayed 
in figure 1 with two or three particles in the intermediate state. We do not have to calculate 
cuts with four particles in the intermediate state since such cuts would run through the charm 
quark bubble and events involving charmed hadrons in the final state are not included on 
the experimental side. 

The reduction of the 2-particle-cut of the first Feynman diagram visualized in figure 1 
to a set of a few so-called master integrals is done by means of the systematic Laporta 
algorithm [23] based on the integration-by-part technique first proposed in [24,25] ^. Since 
the reduction procedure applied to the present problem has already been described in great 

^For an automatized implementation of this algorithm see, e.g., [26]. 
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Figure 1: 2- and 3-particle-cuts of the irreducible b quark selfenergy diagrams with a massive 
charm quark bubble contributing to the b ^ S'j (blue dashed line) and b — sjg (orange 
dashed lines) transitions at 0{a1). Thick lines denote b quarks, thin lines s quarks, wiggly 
lines photons and curly lines gluons. Left-right reflected diagrams are not shown. 



detail in [10, 11], we refrain from repeating it here once more. We just remark that after the 
reduction there remain only three two-loop integrals which do not factorize into a product 
of one-loop integrals. To solve these integrals, we first introduce Feynman parameters in 
the standard way and perform the loop-integrations. At this level, the Feynman parameter 
integrals contain denominators of the form 



{mlPi 



(2.3) 



where Pi and P2 are polynomials in the Feynman parameters. Next, applying the Mellin- 
Barnes representation [27] 

1 f ds oc^ 

-s)V{a + s) , x = mlPi, y = mlP2 , (2.4) 



{x + yY V{a) J 2m 
c 

where the integration contour C runs from —ioo to +ioo such that it separates the poles 
generated by the two F functions, proves very useful because then the integration over the 
Feynman parameters becomes trivial. (For explicite examples of this method, see e.g., [28]). 
Finally, we close the integration contour C by a half-circle with infinite radius at either side 
and sum up the enclosed residues. In this way we managed to obtain solutions for the three 
non-trivial two-loop integrals valid for arbitrary values of m^. 

The remaining 2- and 3-particle-cuts which all contain cut-diagrams with external wave 
function corrections are taken into account by proper insertions of s quark and gluon wave 
function renormalization constants into the corresponding 2- and 3-particle-cuts of irre- 
ducible two-loop b quark self-energy diagrams. A collection of the relevant renormalization 
constants necessary to render our amplitudes finite can be found in appendix A. 

Having briefly described the technical details of our calculation, we turn our attention 
to the results for the function G-j-j{Eq, 11) which we write as an integral over the (rescaled) 
photon energy spectrum: 

G,,{E,,,)= fdz^^Il^, .= ^, .0 = ^. (2.5) 
Jzo dz nib rub 
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In NNLO approximation the photon energy spectrum can be decomposed as follows,'^ 



dz ^ ' ' 47r 

where 



:2.6) 



(2.7) 



Here, A^^^, and A'^v' denote the number of hght {rriq = 0), heavy {rrig = rrib), and virtual 



TTic) quark flavors, that is the total number of quark flavors is Np = Nl + Nh + N] 



V- 



Furthermore, as{fi) is the running coupling constant in the MS scheme, and the numerical 
values of the color factors are given by Cp = 4/3, Ca = 3, and Tr = 1/2. 

The functions if*-^-* and if^^'*) (z = a, na, NL, NH) appearing in (2.6) and (2.7) receive 
contributions from the b — > 57, b — > S'jg, b s^ygg and b S'yqq (g G {u, d, s}, rriq = 0) 
transitions. They can be found in [11]. On the other hand, the function jj^^'^^) is completely 
new. It originates from the 6 — > S7 and b — >■ s'yg transitions with massive charm quark 
bubbles, see figure 1. Our result for this function reads 



124 
"27 



TT 



32 



(7 + TT^) L, 



16 



L.+ 



^-^^^)lnp + /(p))5(l-.) 



16 

+ y (2L^-lnp) 



ln(l - z) 



28 



(2L^-lnp) 



1-z 



--{7 + z-2z^-2{l + z) ln(l - z)) 

+ ^ (7 + z - 2 Inp - ^ (1 + 2) ln(l - ;z) Inp ^ 



(2.^ 



where 
fip) 



^ (162 v/;^ + 70 p3/2 - 36 p2) + ^ p In p + ^ (25 + 27 p^) In^ + ^ ln=^ P 



9 



9 



9 



9 



4 4 

- - (31 + 27p^) ln(l -p)lnp- - (81 + 35p) artanh(7p) In p 
9 9 

- ^ (62 + 81Vp + 35p=^/2 + 54p2) Li2(p) + ^lnpLi2(p) 

8 , 16, , , 5578 172 

+ g (81 + 35 p) Li2(v^) - — Li3(p) + + P • 



(2.9) 



Here, 



m 



p = — I , L^ = \n[ — 



m 



P 



(2.10) 



^In the notation of [11] we have = 4i/(i) and ^ 16i/(2,») for i = a,na,NL,NH. 



5 



Li^iz) = dxLi2(x)/a;, and [•••]+ are plus-distributions defined in the standard way. Note 
that our resuh (2.8) holds for p G (0, oo). The individual 2- and 3-particle-cuts contributing 
to can be found in appendix B. 

In the remainder of this section we comment on the numerical importance of the new 
contribution given in (2.8). To this end we consider the function 



Grr(0,/i) = l + ^ ---^ --L, 



{ (3.55556iVj. - 44.4446) + (21.6168iVj. - 334.803) 

+ 37.8172Nl + h{p)Nv - 2. 16077 A^^ - 519.25o} , (2.11) 

which follows from (2.5) when setting Eq = and performing the integration over z. The 
function h{p) which incorporates the mc-dependence is given by 

^, , 248 2 /1972 16 2\ , 2 ^, , 

Mp) = ^-^+(^-^-^jlnp+3/(p). (2.12) 

Equation (2.11) generalizes our result presented in [10] to include Ny massive charm quarks. 

Denoting the coefficient of (a<i(/i)/(47r))^ in (2.11) by X2{Nh, Nl, Ny), we find (for p = 
nib and mc/nib = 0.26) 

X2(l, 4, 0) = -370.142 , X2(l, 3, 1) = -386.638 , X2(2, 3, 0) = -410.120 , (2.13) 

that is the result with the physical charm quark mass is almost equal to the average of the two 
approximations. From these numbers we conclude that the branching ratio Bt{B —>■ Xs'~f) 
will stay within the errors quoted in (1.2) when implementing the exact mc-dependence of 
the (O7, 07)-contribution. Therefore we do not repeat the interpolation performed in [13] 
which would be necessary to properly implement our new result since it also estimates the 
mc-dependence of the (O7, 07)-contribution. We postpone this until more progress towards 
an improved evaluation of the branching ratio at NNLO has been achieved. 

Before closing this section we mention that a normalization with m^ rather than m1[p)m\ 
is sometimes used in the definition of (2.1). The relation between the MS mass frib{p) and 
the on-shell mass rrtb (including the exact dependence on p) necessary to convert between 
both normalizations can be found in appendix A. 
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A Renormalization constants 



The renormalization schemes apphed in this work are exactly the same as used in [10], and 
the relevant renormalization constants for Ny = can be found in appendix A of that 
reference. In order to obtain the renormalization constants for A'^v' 7^ form those given 
in [10] one has to proceed as follows: (i) set Np, the total number of quarks appearing there, 
equal to Nh + Nl + Nv] (ii) add additional contributions SZ^^, and SZ^^ to the gluon, 
s quark and b quark wave function renormalization constants, respectively. These additional 
contributions are proportional to Ny, and their explicit expressions read 



5Z^^ = CFTRNy}^^ (l + 8L^-4 1np)+ (^y-161np^ L^ + 2ALl 

+ ^ + 28p + ^(5- ISV;^- 30p3/2 + 12 p2) + ^ (2 + 3p) Inp 

+ 2 (2 + 3p^) lnV-4 (1 + 3p^) ln(l - p) Inp 

- 4^ (3 + 5 p) artanh(yp) In p + 8^^ (3 + 5 p) U2{^) 

- 2 (2 + 37^ + 5 p^/^ + 6 p') U,{p) I ' + 0{al) , (A.l) 

with 7 ^ 0.5772 being the Euler-Mascheroni constant and d = 4 — 2e the space-time dimen- 
sion. While the first two expressions in (A.l) can be obtained in a straightforward way, the 
calculation of the analytic expression for 52^^^ is more involved. We checked numerically 
that our result is in agreement with the one given in [29]. We remark that the results in 
(A.l) are independent of the gauge parameter appearing in the gluon propagator. 

Finally, we give the relation between the pole- mass and the MS mass fnb{fi) up to 
two-loops including the contribution of the Ny quarks, 

^ = 1-C.(4 + 6L,)^ 
nif, An 

+ CfIcf (1 + Svr^ In 2 - Stt^ - 12 Cs + 21 + 18 L; 

1111 . 21 o 4 2 185 

h47rMn2 - -TT^ - 6C3 H , ^^^„ 

24 3 2 M I. 



-Ca(^+ 47r2ln2 - -n^ - Q + — + 22 

_ ,143 8 2 52 ^ „ 
+ TrNh { — - -n^ + - L, + 8Ll 



6 3 3 
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/52 71 4 

+ TnNvi - + 8LI + - + 12 p+-7r^l- 3^- 3 p^/^ + p') 

+ 4plnp + 2p^ In^p - 4 (l + p^) ln(l - p) Inp 

- 8i/p(l + p) artanh(yp) Inp + 16^(1 + p) 112(1//)) 

- 4 (1 + v/p + p^/^ + p') Li2(p)) I ' + Oial) , (A.2) 

where is the Riemann zeta-function. We recalculated the term proportional to Ny and 
our finding agrees with [30] (bearing in mind that a factor of 4/3 is missing in front of the 
function A(Mj/M) appearing in equation (17) of that reference, as also observed in [31]). 
The remaining terms in (A.2) have been taken from [30]. 



B Individual cut contributions to //(^'^^) 

The function H^'^'^^'^ introduced in (2.7) receives contributions from the 2-particle-cut of 
the first diagram given in figure 1, as well as contributions from 2- and 3-particle-cuts 
where at least one renormalization constant proportional to A'^v' is present. Denoting these 
contributions by ^^2,Nv,baro)^ ^^2,Nv,ct)^ ^^2,Nv,ct)^ respectively, we have 

^(^'-)(.,/i) = ^^'^-^)(^,/x) + ^f'^^'-)(.,/i) +^(^'^^-^)(^,/z) , (B.l) 
with the individual contributions given by 

188 2 
~9~ ~ 9 



7r2 + [40-161np] + 48 - ^ Inp + ^ InV^ ^ 

7612 80 2/^73 ^ 20 3/2\ „^ 

H p + TT^ 12./P H p^'^ - 8 C3 

81 9 V27 ^ Q J 

(3T6 4 \ 
— - -TT^ -321np + 81n2p] + [120 - 48 Inp] LI + 9QLI 



10 o 8 . . A , . 56 , 2 



+ [j 27 t'^"^^'^]J lnp+— In^p-— Inpln(l-p) 
- ^ (27 - 5 p) artanh(v/p) Inp - ^ (22 + 27y^ - 5 p^/^) Li2(p) 

+ ^ lnpLi2(p) + (27- 5p)Li2(v/p) - y Li3(p) } 5(1 - z) , 

(B.2) 
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- ^ + 28p + vr^ Q - 6Vp - 10 p^/^ + 4p^) + ^ Cs 

/520 16 2 o., ^ r 64 . 256 [2 „ \ 

+ (4 + 6p2) lnV-4 (1 + 3p2) Inp ln(l - p) 

- 4^ (3 + 5 p) artanh(^) lnp-2(2 + 3Vp + 5 p^/^ + 6 p") Li2(p) 



+ 8v/;^(3 + 5p) Li2(v/;^) U(l-z), 



(B.3) 



(2,NV,ct) 



2,P) 



8 1 

-2L^ + lnp)-- 



2 /^56 , \ , 80 . 28, 4, 2 \ 1 



TT 

9 



608 3 



^ - f vr^ - 56 In p + 8 In - - 48 In p ) 



f - 2 vr^) In p - ^ InV + ^ In^ P + ^ Cs M( 1 - ^) 



16 

+ y (2L^-lnp) 



8 



ln(l - z) 
l-z 



28 



(2L^-lnp) 



l-z 



-{7 + z-2z^ -2(1 + z) ln(l - ^)) 

4 8 
+ - (7 + ^-2^2) Inp - - (1 + 2)ln(l - 2)lnp. 
3 3 



(B.4) 



The three contributions given above are by themselves independent of the gauge parameter 
entering the gluon propagator. 
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